Patterned Deposition of Particles in Spatio-temporally Driven Lattices 
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We present and analyze mechanisms for the patterned deposition of particles in a spatio- 
temporally driven lattice. The working principle is based on the breaking of the spatio-temporal 
translation symmetry, which is responsible for the equivalence of all lattice sites, by applying mod- 
ulated phase shifts to the lattice sites. The patterned trapping of the particles occurs in confined 
chaotic seas, created via the ramping of the height of the lattice potential. Complex density profiles 
on the length scale of the complete lattice can be obtained by a quasi-continuous, spatial deformation 
of the chaotic sea in a frequency modulated lattice. 
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Introduction The desire of understanding the fascinat- 
ing mechanism that allows living organisms to extract 
mechanical energy from thermal fluctuations evoked elab- 
orate theoretical and experimental investigations [1]. A 
directed flow of particles can in principle arise from ther- 
mal fluctuations, if an appropriate symmetry-breaking 
and suitable time-correlations are present [2[. The spatio- 
temporal symmetries to be broken lead to the existence 
of a counter-propagating partner to each trajectory [S^. 
Ratchet transport has consequently been shown to occur 
in the nonequilibrium dynamics of driven Hamiltonian 
systems with mixed phase space . A first experimen- 
tal realization of an electron ratchet in a semiconductor 
nanostructure [7] showed the invertibility of the direction 
of the current via the temperature. The full range of con- 
trollability of the ratchet-transport became clear in ex- 
periments in Josephson junction- arrays [8] and especially 
in the remarkable realizations of ratchets with atoms in 
time-dependent optical lattices |9I-I1 1\ . On the theoretical 
side, most work has focused on the dynamics of a single 
barrier potential [12] and on setups where the spatial and 
temporal dependencies of the underlying forces decouple 
B m, [Tsl, 13 • is o^ly recently that a driven lattice with 
a spatially varying driving law has been theoretically ex- 
plored [15I] and the feasibility of a local engineering of 
phase space has been demonstrated. Motivated by this 
enhanced control achieved in phase modulated lattices 
[15], we address here the question: Is it possible to ex- 
ploit spatio-temporally varying driving laws in order to 
realize a mechanism for selective trapping of particles in 
lattices thereby leading to certain patterns of their spa- 
tial deposition? 

In the present work, we demonstrate and analyze a mech- 
anism that provides us with specific pattern-like distri- 
butions of particles in a lattice of periodically driven 
potential barriers, starting from a homogeneous initial 
distribution. Considering a lattice with a certain bro- 
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FIG. 1. (Color online) Schematic illustration of the driven 
lattice and the trapping mechanism: Before ramping the po- 
tential (a), after ramping (b) and at the end of the simulation 
(c). Blue particles are trapped, red ones are diffusive. 



ken spatio-temporal symmetry and ramping (i.e. instan- 
taneously increasing) the lattice potential after a given 
time, we observe that particles are trapped on specific 
sites of the unit cell, while the other sites are diffusively 
emptied. Here, the periodicity of the resulting particle 
distribution follows the length scale of the unit cell. We 
also present a second mechanism for the creation of so- 
phisticated particle distributions following certain glob- 
ally varying spatial patterns. 

Setup We consider a system of noninteracting classical 
particles in a finite one-dimensional lattice of laterally 
oscillating potentials with equal height described by 
the following single particle Hamiltonian: 
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Therein, m, x and p denote the mass, position and mo- 
mentum of the particle and a, a;, f{t) = cos(cjt + (pi) 
represent the amplitude, frequency, width and driving 
law of the oscillating barriers. L is the equilibrium dis- 
tance between two barriers and is the phase of the 
z-th barrier. A sketch of the lattice is shown in Fig. la. 
The space between the z-th and the {i + l)-st barrier 
will be referred to as the i-th lattice site. As a parti- 
cle reaches the ±(A/'/2 + l)st barrier it is considered to 
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be lost. Obviously the system is unbiased, i.e. the time 
averaged net force vanishes. With the scaling transfor- 
mations t ^ t' \= uot] X ^ x' \= x/a and additionally 
H'{x' ^t') := the number of parameters can be re- 
duced by three. Therefore, we may choose m = I = 1 
and L = 2 + 2a + / without loss of generality. 

Spatio-temporal symmetry breaking What kind of 
spatio-temporal driving is necessary in order to achieve 
selective trapping into certain sites of a lattice? Obvi- 
ously, our aim cannot be achieved with a 'global' (i.e. site 
independent) driving law, since then trajectories would 
obey the same dynamics on every lattice site. In Ref. 
[Hj local changes in phase space were obtained in a lat- 
tice by employing linear phase shifts = i • (j). They 
are however not sufficient to lead to selective trapping of 
particles, since the following symmetry still holds: 



H{x ^L,t- = H{x, t) 



(2) 



In lattices satisfying the symmetry (|2j) a trajectory which 
is at time t at a position x has the same future as a tra- 
jectory which is at t — (j)/uj at the position x + L. That 
is, if particles are trapped on site z, those on site z + 1 
are also trapped. For this reason, breaking the symme- 
try (j2j) of the Hamiltonian is a necessary requirement for 
selective trapping. 

Following the above arguments we choose a system with 
a unit cell consisting of three barriers, i.e. 0^+3 = (j)i and 
nonlinear phase shifts (^o, (/>2) = (0, tt/IO, Stt/IO). 
Initially, our lattice consisting of = 4000 barriers, pos- 
sesses a low potential height V = 1 and we choose a uni- 
formly distributed ensemble of Np = 10^ identical, non- 
interacting classical particles with randomly distributed, 
small velocities (corresponding to initial conditions in the 
chaotic sea of the underlying phase space). The parti- 
cle dynamics is followed until a system-specific long-time 
transient and a quasi-stationary distribution in the lat- 
tice are reached. Then at t^ := 2 x 10^ (which corresponds 
for each particle approximately to the same number of 
collisions with a barrier), the potential is ramped (in- 
stantaneously increased) to the larger value V = 9 and 
the dynamics is followed until tf = 10^. 

Results Fig. [2^ shows an extract of the time-evolution 
of the particle density in the lattice. Initially, the particle 
distribution is uniform. After ramping the potential, we 
observe that the particle number remains approximately 
constant on every third lattice site, while the other sites 
are emptied (see Fig. [TJd). Asymptotically, the distri- 
bution becomes stationary (see Fig. [It) and only every 
third lattice site is occupied by particles whereas the in- 
between sites are empty. Analogously, Fig. [SJd reveals, 
that ramping the potential from V = ltoF = 15 results 
in trapping on two of three sites of the unit cell. 
We now turn to a quantitative discussion of the under- 
lying process and address the questions of efficiency and 
predictability of the selective trapping mechanism. Tab.|I] 
shows the final particle distribution at t = tf ^ and the 
first three columns list the fraction of particles (in per- 
cent), which are trapped on the first /second/third site 
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FIG. 2. (Color online) Representative extract of the time- 
evolution of the density of particles: The particle number 
(color) per site (density) is plotted versus time (logarithmic 
scale) and site number. The potential is ramped at tr = 2 x 
10^ from y = 1 to y = 9 (a) and to y = 15 (b). (a = cj = 1). 



of the unit cell, taken with respect to the total number 
of particles which are in the lattice at t = tr. We note, 
that trapping on the first site of the unit cell corresponds 
of course to trapping on the sites i = .., — 2, 1, 4, 7, .. in 
the extended lattice. Without ramping the potential, we 
observe from the V = 1-row, that at t = tf the lattice 
is completely emptied. Ramping from ]/ = 1 to F = 9 
(to V = Ib/V = 24), we see that the fraction of particles 
trapped on the first site (the first two sites/all three sites) 
of the unit cell exceeds the fraction of particles trapped 
on the other (undesired) sites by about a factor 100. The 
small numbers of particles trapped on the undesired sites 
are a consequence of the impossibility to distinguish long 
transient trapping from permanent trapping. The sec- 
ond last column of Tab. |I] displays the overall fraction 
of particles which is trapped at tf. This number con- 
verges for high barrier potentials to the system specific 
limit of about 82% of all available particles, shown in 
the case V = 100. Finally Tab. |I] shows the number of 
particles, which are neither trapped nor have left the lat- 
tice at the end of the simulation. Since this fraction of 
particles which are asymptotically diffusive is small, it is 
confirmed that the obtained pattern-like deposition be- 
comes indeed stationary. Note that all the relative num- 
bers provided in Tab. U are independent of the specific 
initial particle distribution. The ratio of the fluctuations 
(standard deviation) and the average number of particles 
per site is expected to decrease asymptotically for large 
particle numbers, according to with a system 

specific proportionality constant. Our simulation of 10^ 
particles yields fluctuations of the particle number per 
site of 2.4% for V = 9. 

In the following, we show that the mechanism of selective 
trapping can be traced back to the Poincare surface of 
section of the system. We firstly discuss how the breaking 
of the symmetry ([2]) leaves its hallmarks in the underlying 
phase space. In Fig. [3^ we observe that for the low barrier 
potential V = 1^ the chaotic sea represents a connected 
area. Therefore, diffusion occurs all-over the lattice. For 
a higher barrier potential V = 9, according to Fig. [3]d a 
part of the chaotic sea is separated (region T) by regular 
spanning curves from the main chaotic sea which corre- 
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spends to diflFusion over the whole lattice (region F). The 
particles in region T never leave the corresponding lat- 
tice site, i.e. they are trapped. Simultaneously, particles 
in region F follow a diffusive motion until they reach the 
border of the lattice and leave it. For a potential height 
F = 15, Fig. [3fc shows two isolated chaotic seas (regions 
Ti, T2). Starting again with V = 1 and ramping to 
F = 15, trapping thus occurs on two sites of each unit 
cell, as we have demonstrated in Fig. [JJd. Finally, for 
V = 24, the chaotic sea, (Fig. [3]i) suggests trapping on 
all lattice-sites. Every lattice site can be regarded as a 
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FIG. 3. (Color online) Stroboscopic Poincare surfaces of sec- 
tion (PSOS) for the unit cell at times t mod 27r = for 
different values of the barrier potential 1/ = 1, 9, 15, 24 for (a) 
to (d), respectively. Components corresponding to a local- 
ized chaotic sea (blue) are denoted by "T", those which are 
extended over the whole lattice by "F" (black). Red curves 
show regular orbits, {a = u = 1). 

Fermi-Ulam- model (FUM) with finite wall poten- 

tial. Accordingly, the trapping mechanism can be under- 
stood on a physically intuitive level: In the FUM, due to 
the formation of correlations, there are regular spanning 
curves, which prohibit unlimited energy growth. If the 
kinetic energy corresponding to the first invariant span- 
ning curve (FISC) lies completely (i.e. for all phases of 
the oscillation and all coordinates) below the barrier po- 
tential height, particles below the barrier will not be able 
to overcome it. In contrast to this, particles in the re- 
gion F will escape from the lattice after some time. In 
the investigated symmetry-broken lattice, the unit cell 
can be considered to consist of three different finite- wall- 
FUMs. Thus, particles will be trapped on those sites of 
the unit cell which have FISCs which are completely be- 
low V. Ramping the potential to V = 9 {V = 15), we 
encounter trapping in the unit cell on the site with the 
lowest (two lowest) FISC(s), since the FISC leads to a 
confined chaotic sea for the corresponding site(s). 
According to Tab. Ueven for a very large potential V the 
number of trapped particles cannot exceed 82%. This is 
due to the fact that particles which are at the location of 
a barrier at the moment when the potential is increased, 
gain the potential energy V{t > t^) — V{t < t^) and 



leave the lattice either on regular orbits or via diffusion 
above the FISC. Since the particle distribution for t < 
is close to uniform, the fraction of maximally trappable 
particles can be estimated as L — l/L. For very small 
barriers / ^ 0, the overall fraction of trappable particles 
is expected to approach 100% which was confirmed by 
our numerical simulations. The typical velocities of the 
FISC on a specific site strongly increases with the phase 
difference of the adjacent barriers for not too large values 
of this difference. We can thus develop a procedure that 
yields a lattice which contains particles arranged accord- 
ing to a desired pattern. We might for example wish as 
a final distribution to have particles only on every sec- 
ond lattice site. To this end, a period-4-lattice with a unit 
cell consisting of two barriers with small and two barriers 
with large phase shifts serves the purpose. For example, 
by choosing (0o, ^1, 02, ^3) = (0, ^, ^, 0) and ramping 
from V = lioV = 100, we obtain a final distribution of 
(0.055/20.90/0.003/20.82; 41.77; 3.98) following the no- 
tation of Tab.m Also the relative occupation of the occu- 
pied sites can be changed by an additional modulation of 
the barrier amplitudes a ^ a^; = a^+a; z = 0, 1, 2;. 
Choosing ai = 2.0; a2 = 0.5; as = 1 and ramping 
from V = 1 to V = 100 in the period-3-case, we obtain 
(lst/2nd/3rd; t.f.; d.f.) = (36.30/18.20/27.96; 82.44; 0.20) 
and for ai = 0.5; as = 2.0; as = 1: (23.83/37.98/20.69; 
82.47; 0.20). This demonstrates the tunability of the rel- 
ative occupation probability of the lattice sites. 
The mechanism of patterned deposition is expected to 
be stable with respect to small thermal fluctuations of 
the ensemble of particles, because the breaking of the 
symmetry (j2j) is not affected by the latter. However 
they will cause a certain transmissibility of the spanning 
curves, because the energy change provided by the ther- 
mal fluctuations enables particles to overcome the span- 
ning curve. As more and more spanning curves occur 
with increasing V ^ this leaking rate can be tuned to an 
arbitrary small value. 

Since the ratchet effect can occur for both, particles 
in spatio-tempor ally-driven potentials, and thermally 
driven particles in a static potential, it is interesting to 
ask whether the mechanism of patterned deposition can 
also be obtained in thermally driven setups. For a sys- 
tem, consisting of two species, an enormous flexibility, 
concerning direction and magnitude of the ratchet cur- 
rent was demonstrated ^18r.20j , thus suggesting that peri- 
odic clustering can also occur in a thermally driven setup. 



Global patterns of particle deposition We explore now the 
question, whether a selective trapping of particles is pos- 
sible, on the length scale of the complete lattice, i.e. fol- 
lowing patterns beyond the size of a few-site unit-cell. 
It turns out, that a frequency- modulated lattice is the 
key to achieve global modulations of the density. Let us 
consider a uniformly distributed ensemble of 10^ parti- 
cles with initially low, randomly distributed velocities in 
a lattice consisting of = 10^ barriers possessing a site 
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TABLE I. Propagation of 10^ particles, for t — 10^ time steps 
in a lattice consisting of 4000 individual barriers (parameters: 
Fig. [21). At the ramping time tr = 2 x 10^, 19.64% of all lO'^ 
particles have left the lattice. First column: Barrier-potential 
V{t > tr). "site 1(2,3)": Conditional probability for a particle 
which was still in the lattice at tr, to be trapped on the first 
(second, third) site of the unit cell at t = tf. "total fraction": 
Fraction of trapped particles, related to the total number of 
particles which are in the lattice at tr- The total fraction of 
diffusive particles in per cent at t = 10^ is also shown. 



V 


site 1 


site 2 


site 3 


total fraction 


diffusive fraction 


1 


0.00 


0.00 


0.00 


0.00 


0.36 


9 


26.54 


0.15 


0.26 


26.95 


0.38 


15 


26.66 


23.38 


0.27 


50.32 


0.37 


24 


26.66 


23.38 


27.67 


77.71 


0.77 


100 


26.66 


23.38 


31.93 


81.97 


0.37 



dependent frequency oo ^ ooi := 1 -\- sin(i/1500). We in- 
tegrate the dynamics for V = 1 up to tr = 2 x 10^ and 
subsequently ramp the barriers to V = 120, followed by 
a further propagation to tf = 10^. 

Fig. m shows the corresponding evolution of the particle 
density at four different times. Initially, the ensemble 
is uniformly distributed. For t <C tr the emergence of 
a (quasi-) periodic global density modulation on large 
length scales is conspicuous. The peaks develop in those 
regions of the lattice where the barrier frequency is lower. 
For t < tr the observed modulation is more pronounced. 
Finally, for t^^ the high-frequency regions of the lat- 
tice are completely emptied. The resulting density dis- 
tribution has become stationary. 
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FIG. 4. (Color online) Snapshots of the time- evolution 
of the particle-density in the lattice with a spatially vary- 
ing frequency uji — 1 -\- sin(i/1500) (Inset: uji — 1.2 + 
0.25sin(Vl500) - 1.6 • \i\/N) (a = 1). 

Let us investigate the mechanism for the time-evolution 
of the particle density leading to a localization of the par- 
ticle distribution in the regions of the lattice with slowly 
oscillating barriers. Here, we consider the behavior of 
the FISC on a specific lattice site. In the high frequency- 
regions of the lattice, the FISC occurs at higher veloci- 
ties than in the low- frequency regions. This leads to two 
important effects: Firstly, we have a relatively high av- 
erage particle velocity in the high frequency-regions, cor- 
responding to a short average dwell time and resulting 
in a low particle density when compared to regions with 



slowly oscillating barriers. Since the barrier frequency 
varies quasi-periodically over the lattice, the resulting 
particle density varies likewise with the same periodic- 
ity. Secondly, the maximal kinetic energy a particle can 
achieve on a lattice site increases with the typical ve- 
locity of the FISC. Thus, after ramping the potential to 
an appropriate value, following the above discussion of 
the selective localization in a few-site unit cell, particles 
will be trapped on those sites, where the FISC is located 
completely below the potential height. The latter occurs 
in the low- frequency regions of the driven lattice. 
On basis of the above-developed understanding, we can 
now further design particle distributions by appropriately 
choosing the barrier frequencies: For instance, the final 
particle distribution for the superposition of a periodic 
frequency- modulation and a piecewise constant gradient, 
declining from the center of the lattice to both edges 
uJi = 1.2 + 0.25sin(z/1500) - 1.6 • \i\/N, is shown in the 
inset of Fig. |4j The density reflects the expectation of a 
particle accumulation in the low- frequency regions. Note 
that for the continuous manipulation of the FISC it is 
necessary to have a small frequency change from barrier 
to barrier. Otherwise, particles in the fast regions could 
move on regular orbits traversing the low frequency re- 
gions, instead of being captured in the chaotic sea in the 
course of the ramp. 

Conclusion We have demonstrated and analyzed mech- 
anisms for the patterned deposition of particles in a lat- 
tice by applying spatially varying driving laws and a 
corresponding ramp of the potential height. The sym- 
metry which leads to the long-term equivalence of all 
lattice sites has been broken by appropriately chosen 
phase shifts of the barriers. This results in trapping in 
a patterned series of confined chaotic seas. Employing a 
frequency modulated lattice we demonstrated that even 
patterned density profiles on the length scale of the com- 
plete lattice can be achieved. A possible experimental 
setup to demonstrate the mechanisms of selective trap- 
ping is provided by cold atoms in optical lattices. Oscil- 
lating mirrors and acousto-optic modulators can be used 
to obtain a laterally oscillating potential landscape. Su- 
perimposing laser beams of different frequency (ampli- 
tude) and adjusting the relative phases by controlling 
the positions of the light sources on a nm-scale we ex- 
pect to arrive at setups which show the demonstrated 
mechanism of patterned deposition. The ramp leading 
to the final confinement in the individual sites can be 
implemented by an instantaneous switch of the ampli- 
tude of the laser (s). An alternative physical system to 
realize our model is provided by multi-layered semicon- 
ductor heterostructures, where applied AC voltages can 
be used to drive each layer individually. Potential ap- 
pHcations include selective loading of a lattice with cold 
atoms, thus permitting local addressing for quantum in- 
formation processing. It might be also of relevance to 
atom lithography. 
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